arXiv: 1508.06284vl [astro-ph.GA] 25 Aug 2015 


Accepted to ApJ on August 14, 2016 

Preprint typeset using style emulateapj v. 5/2/11 


MICROWAVE CONTINUUM EMISSION AND DENSE GAS TRACERS IN NGC 3627: COMBINING JANSKY 

VLA AND ALMA OBSERVATIONS 

Eric J. Murphy\ Dillon Dong^, Adam K. Leroy^, Emmanuel Momjian'‘, James J. Condon®, George Helou®, 
David S. Meier^’'‘, Jurgen OtU, Eva Schinnerer®, Jean L. Turner® 

Accepted to ApJ on August I 4 , 2016 

ABSTRACT 

We present Karl G. Jansky Very Large Array (VLA) Ka band (33 GHz) and Atacama Large Millime¬ 
ter Array (ALMA) Band 3 (94.5 GHz) continuum images covering the nucleus and two extranuclear 
star-forming regions within the nearby galaxy NGG3627 (M66), observed as part of the Star For¬ 
mation in Radio Survey (SFRS). Both images achieve an angular resolution of <2", allowing us to 
map the radio spectral indices and estimate thermal radio fractions at a linear resolution of <90 pc 
at the distance of NGC 3627. The thermal fraction at 33 GHz reaches unity at and around the peaks 
of each Hii region; we additionally observed the spectral index between 33 and 94.5 GHz to become 
both increasingly negative and positive away from the peaks of the Hii regions, indicating an increase 
of non-thermal extended emission from diffusing cosmic-ray electrons and the possible presence of 
cold dust, respectively. While the ALMA observations were optimized for collecting continuum data, 
they also detected line emission from the J = 1 ^ 0 transitions of HGN and HCO'*'. The peaks of 
dense molecular gas traced by these two spectral lines are spatially offset from the peaks of the 33 and 
94.5 GHz continuum emission for the case of the extranuclear star-forming regions, indicating that our 
data reach an angular resolution at which one can spatially distinguish sites of recent star formation 
from the sites of future star formation. Finally, we find trends of decreasing dense gas fraction and 
velocity dispersion with increasing star formation efficiency among the three regions observed, indi¬ 
cating that the dynamical state of the dense gas, rather than its abundance, plays a more significant 
role in the star formation process. 

Subject headings: galaxies: individual (NGC 3627) - galaxies: star formation - Hii regions - ISM: 
molecules - radio continuum: general 


1. INTRODUCTION 

Stars form out of small pockets of dense molecular 
gas, with massive stars that are unambiguously young 
often embedded behind thick layers of dust. The dense 
star-forming gas is selectively traced by transitions with 
high critical densities, like the J = 1 —>■ 0 transitions of 
HCN or HCO< (i.e., ricrit ^ lO^crn”^). Likwise, recently 
formed stars are best studied by tracers such as free- 
free continuum at A < 1 cm or hydrogen recombination¬ 
line tracers that are robust against extinction and selec¬ 
tively sensitive to the massive young stars that live only 
a short time after their birth. To constrain the complete 
physical process of star formation, ideally both tracers 
would be observed on scales matched to individual star- 
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forming regions. However, because of the faintness of 
both types of emission and the need for wide frequency 
coverage to isolate free-free emission from contaminants, 
this sort of study has been challenging before the cur¬ 
rent generation of radio and millimeter-wave telescopes, 
especially the Karl G. Jansky Very Large Array (VLA) 
and the Atacama Large Millimeter/submillimeter Array 
(ALMA). In this paper, we combined new observations 
from ALMA and the VLA to report one of the first extra- 
galactic cloud-scale comparisons of dense gas, traced by 
the high effective density transitions HCN (J = 1 —0) 
and HCO< (J = I —i 0), and recent star formation 
traced by A < I cm free-free emission. 

This paper represents the first results of a larger 
project to combine ALMA and the VLA to study contin¬ 
uum emission from nearby galaxies at frequencies rang¬ 
ing between 3 — 100 GHz. Emission at frequencies span¬ 
ning 30 — 100 GHz is expected to be dominated by free- 
free emission from Hii regions, providing a highly robust 
measure of massive star formation activity unbiased by 
dust. Unfortunately, this emission component is energet¬ 
ically weak, making detections difficult and time consum¬ 
ing even for bright objects. To date, observations in this 


frequency ran ge have been largely limited to Galactic Hii 
regions (e.g., iMezger fc Henderson] |196 7ll. ii earbv dwarf 


irregular galaxies fe.g.. Klein & GraevellOSOh. galaxvnii- 

clei ('e.e.. iTurner & Hd 

19831 1994). nearbv starbursts 

le.g.. iKlein et al.l 119881 

Turner & Hd 1198511. and snoer 

star clusters within nearby blue compact dwarfs (e.g., 

ITurner et al.l 119981: iKobnlnickv & JohnsonI 1199911. How- 

ever, with ALMA now 

online, 3 mm bolometer arrays 
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Table 1 

Source Names and Imaging Summary 


Source 

R.A. 

(J2000) 

Decl. 

(J2000) 

Beam (33 GHz) 

C733 GHz 
(/xJybm-1) 

Beam (94.5 GHz) 

^^94.5 GHz 
(/iJybm-1) 

NGC 3627 

11 20 15.0 

-1-12 59 30 

2'.'83 X T.'83 

26 

T.'89 X 1'.'66 

32 

NGC 3627 Enuc. 1 

11 20 16.2 

+11 20 16 

2'.'45 X 2'.'03 

20 

T.'88 X 1'.'67 

32 

NGC 3627 Enuc. 2 

11 20 16.3 

+11 20 16 

2'.'55 X 2'.'08 

25 

1'.'89 X 1'.'67 

32 


such as MUSTANG (iDicker et alJl200^ . along with re¬ 
cent improvements to the backends of existing radio tele¬ 
scopes, such as the Caltech Continuum Backend (CCB) 
on the Robert C. Byrd Green Bank Telescope (GBT) 
and the Wideband Interferometric Digital ARchitecture 
(WIDAR) correlator on the VLA, the availability of in¬ 
creased bandwidth is making it possible to conduct in¬ 
vestigations for larger sampl es of extragalytic objects 
at frequencies >30 GH z fe.s.. iMurphv et 311120111 [20121 : 
iNikolic fc Boltonll2012ll . 

Taking advantage of this new ca pability, the Star For - 
mation in Radio Survey (SFRS; see lMurphv et alJl2012ll . 
targets 118 star-forming regions (56 nuclear and 62 ex- 
tranuclear) in 56 nearby galaxies (d < SOMpc) that 
had been observed at infrar ed (IR) and optical w ave- 
lengths as p art of the SINGS (iKennicutt et al.ll2003ll and 
KINGFISH (jKennicutt et alji2nilh legacy programs. Of 
these, 112 (50 nuclei and 62 extranuclear regions) have 
6 > —35° and are thus observable with the VLA, while 
54 are observable with ALMA (i.e., S < 30°). We ob¬ 
served 9 of these regions at frequencies near 95 GHz using 
ALMA’s Band 3 receiver during the second early science 
campaign (“Cycle 1”). All 9 of these targets had already 
been observed by the VLA. 

Thanks to ALMA’s strong multiplexing capabilities 
and excellent sensitivity, we were also able to make a si¬ 
multaneous search for bright line emission near 95 GHz. 
Specifically, the bandpass always covers the rest-frame 
J = 1 —>■ 0 transition of HCO'*' and often includes 
the J = 1 —>■ 0 transition of HCN. These are two of 
the standard extragalactic tracers of dense molecu lar gas 
(|Gao fc Solomonir2004lGracia-Carpio et al.ir2006ll . Thus 
our combined VLA and ALMA observations yield an un¬ 
biased view of the ionizing radiation emitted by massive 
young stars and a snapshot of the kind of bright, dense 
structure that might give birth to these regions. 

In this paper, we report on the combination of VLA 
and ALMA data for three regions in the nearby galaxy 
NGC3627 (M66), having a kinematic local standard 
of rest velocity i^lshk = 726 ± 3km s“^. Along with 
NGC3628 and NGC3623, NGC3627 makes up the well- 
known Leo Triplet galaxy groim. It is classified with 
a morphological ty pe of SABtfl, located at a distance 
of 9.38 ± 0.35 Mpc (iFreedman ^ al.l 1200111 , and hosts a 
Seyfert2 AGN (iHo et al.lll997t iMoustakas et al.l 1201011 . 
An Hi plume extending R:i50' indicates that the two 
largest spirals in the grou p, NGC 3627 and 3 628, have in¬ 
teracted in the past ie.g.. lHavnes et al1ll979[ l. NGC 3627 
has comparable amounts of atomic and molecular gas 
(jHelfer et al.ll200,3l : IWalter et al.ll200i^ , which is a rather 
high molecular gas fraction compared to other local star- 

^ Kinematic LSR velocity and morphological type 
taken from the NASA Extragalactic Database (NED; 
http: / / nedwww.ipac.caltech.edu). 


l’^2.7kpc 



Figure 1. The locations of the VLA (dashed line) and ALMA 
(solid line) pointings for each region in NGC 3627 overlaid on a 
Spitzer A = 8 pm grayscale image. The circles outline the FWHMs 
of the VLA and ALMA primary beams at 33 and 94.5 GHz, re¬ 
spectively. 

forming galaxies ie.g.. lSaintonge et al.ll2011h . It has been 
suggested that the high H 2 /H 1 mass ratio in NGC 3627 is 
the result of the tidal interaction with NGC 3628, since 
this galaxy has stripped mu ch of the Hi originally in 
NGC 3627 (iZhang et al.lUb^ . 

The paper is organized as follow: In §2 we describe the 
data as well as the analysis procedures used in the present 
study. In §3 we present our results and discuss their 
implications. Our main conclusions are then summarized 
in §4. 

2. DATA AND ANALYSIS 

We targeted the nucleus and two extranuclear star¬ 
forming regions in NGC 3627 with both ALMA and the 
VLA (Table [J), achieying a common resolution of , 
which is ~90pc at the distance of NGC 3627. In Fig¬ 
ure [T] we identify the locations of each region targeted 
in NGC 3627, showing the full width at half maximum 
(FWHM) of the VLA and ALMA primary beams at 33 
and 94.5 GHz on a Spitzer A = 8 pm greyscale image. 
The fields were selected to target individual massive star¬ 
forming regions, and the A = 8 pm image shows that each 
field is centered on a substantial concentration of bright 
polycyclic aromatic hydrocarbon (PAH) emission, indi¬ 
cating a combination of intense radiation and abundant 
small dust grains. For the case of NGC 3627, this yields a 
sample that includes the nuclear disk, an isolated Hii re¬ 
gion (Enuc. 1), and the star-forming complex at the end 
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Figure 2. Continuum images of the three regions at 33 and 94.5 GHz are shown in the first and second rows, respectively. Coordinates are 
in J2000. The restoring beam size (FWHM) is illustrated in the upper right corner of each panel. At the distance of NGC 3627, 1" linearly 
projects to ss45.5pc on the sky. CO contours are overlaid on each greyscale image as black dashed lines. The resolution of the CO images 
is far lower than in our ALMA or VLA images, being 7'.'3 X 578. HCN contours (cyan) are overlaid on the 33 GHz grayscale images in 
the top three panels, while HCO+ contours (red) are overlaid on the 94.5 GHz grayscale images in the bottom three panels. In all cases, 
contour lines are drawn using a linear scaling starting from the 2cf value. As discussed in the text, the full velocity range of HGN emission 
towards the nucleus is cutoff by the bandpass at ulsrk ~ 800kms“^, and therefore not shown. 


of a bar (Enuc. 2), allowing us to compare star forma¬ 
tion and molecular gas properties for three very different 
environmental conditions. 

We worked with four data sets: VLA continuum ob¬ 
servations at Ka band (29-37GHz), ALMA continuum 
observations centered at w95 GHz, ALMA line observa¬ 
tions of the J = 1 —>■ 0 transition of HGO’*' and HCN, 
and archival BIMA Survey of Nearby Galaxies (SONG 
iHelfer et al.ll200^ line observations of the J = 1 —>■ 0 
transition of CO. 

2.1. VLA Ka-Band Data 

Details on all of the SFRS VLA Ka-band survey obser¬ 
vations and data reduction can be found in E.J. Murphy 
et al. (2015, in preparation). Here, we present data 
only for NGC 3627. D-configuration observations were 
obtained in November 2011 (VLA/llB-032) and March 
2013 (VLA/13A-129). For the first round of observa¬ 
tions, the 8-bit samplers were available, yielding 2 GHz 
of simultaneous bandwidth, which we used to center 
1 GHz wide sub-bands at 32.5 and 33.5 GHz. For the 
latter run, the 3-bit samplers became available, yielding 
8 GHz of instantaneous bandwidth in 2 GHz wide sub¬ 
bands centered at 30, 32, 34, and 36 GHz. In both cases, 
3C286 was used as the flux density and bandpass cali¬ 
brator, while J 1118-1-1234 was used as the complex gain 
and telescope pointing calibrator. To reduce the VLA 
data, we used the Common Astronomy Software Appli¬ 


cations fCASA: [McMullin et al.l[200^ and followed stan¬ 
dard procedures. 

2.2. ALMA Band-3 Data 

In December 2013 as part of ALMA’s Cycle 1 observ¬ 
ing campaign, we obtained data for a subset of 9 SFRS 
sources chosen to be sufficiently bright at 24 /im and close 
enough on the sky that multiple sources could be ob¬ 
served in a single scheduling block. The array configura¬ 
tion was chosen to match the <2"synthesized beam of the 
VLA 33 GHz data. While the goal of the ALMA observ¬ 
ing program is to deliver continuum images at r:!95 GHz, 
we set the local oscillator frequency to 94.5 GHz, center¬ 
ing the four 1.875 GHz wide spectral windows at 87.5, 
89.5, 99.5, and 100.5 GHz, to cover the rest frequencies 
of HCN (J = 1 ^ 0)/88.6318GHz, HCO+ (J = 1 
0)/89.1885GHz, HNC (J = 1 ^ 0)/90.6636GHz, and 
H40q!/99.0229 GHz. Like the VLA data, the ALMA data 
were reduced and calibrated using CASA following stan¬ 
dard procedures as part of the ALMA quality assurance 
process. Here we present the nucleus and two extranu- 
clear star-forming region in NGC 3627. The remaining 
ALMA sources, along with corresponding VLA data, will 
be presented together in a forthcoming paper. 

2.3. Archival CO (J = 1 ^ 0) Data 

To investigate how the amount of dense star-forming 
gas compares to the total molecular gas reservoir, we 
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Table 2 

VLA and ALMA Continuum Properties 


Source 

R.A. 

(J2000) 

Decl. 

(J2000) 

Sp 

(mjybm 

Si 

(mjy) 

dyi X 6m 

Tb 

(K) 

33 GHz Components 

NGC 3627 

NGC 3627 Enuc. 1 
NGC 3627 Enuc. 2 

11 20 15.01 
11 20 16.30 
11 20 16.46 

-1-12 59 29.5 
-1-12 57 49.1 
-1-12 58 43.5 

1.70 ±0.06 
0.77 ±0.03 
0.98 ±0.04 

2.09 ±0.08 
1.66 ±0.09 
5.35 ±0.27 

1'.'21 ± 0'.T5 X 0797 ± O'.TO 
2'.'52 ± 0'.'23 X 2725 ± 0727 
6795 ± 0730 x 3712 ± 0717 

1.98 ± 0.34 
0.33 ± 0.05 
0.28 ±0.02 

94.5 GHz Components 

NGC 3627 

NGC 3627 Enuc. 1 
NGC 3627 Enuc. 2 

11 20 15.02 
11 20 16.29 
11 20 16.44 

-1-12 59 29.6 
-1-12 57 49.2 
-1-12 58 43.5 

1.22 ±0.07 
0.63 ±0.05 
0.68 ±0.05 

1.58 ±0.11 
1.13 ±0.11 
3.23 ±0.25 

I7l0 ± 0716 X 0777 ± 0729 
1795 ± 0724 X 1724 ± 0722 
5766 ± 0736 x 1785 ± 0715 

0.252 ± 0.102 
0.063 ± 0.015 
0.042 ± 0.005 


make use of J = 1 - A 0 CO data take n as part of the 
BIMA SONG survey (|Helfer et al.ll200^ . The rms noise 
of the CO channel map is 41mJybm“^ in a 10kms“^ 
channel. The resolution of the CO map is significantly 
coarser than our ALMA and VLA data, having a syn¬ 
thesized beam of 7"3 x 5'.'8. 

2.4. Continuum Imaging and Photometry 

Calibrated VLA and ALMA measurement sets for each 
source were imaged using the task CLEAN in CASA. The 
Ka-band images contain data from both sets of observa¬ 
tions, but are heavily weighted by the 13A semester ob¬ 
servations as those include significantly more data. The 
mode of CLEAN was set to multifrequency syn thesis (mfs; 
iConwav et ahlfTO^ iSault &: Wieringal 1199411 . We chose 
to use Briggs weighting with robust=0.5, and set the 
variable NTERMS=2, which allows the cleaning procedure 
to also model the spectral index variations on the sky. To 
help deconvolve extended low-intensity e mission, we took 
advantage of the mult iscale clean option (|Cornwellll2008t 
iRau fc Cornwelll[2011ll in CASA, searching for structures 
with scales and 3 times the FWHM of the synthe¬ 
sized beam. A primary beam correction was applied us¬ 
ing the CASA task impbcor before analyzing the im¬ 
ages. The primary-beam-corrected continuum images at 
33 and 94.5 GHz for each of the three targeted sources 
are shown in Figure [5] The synthesized beamwidths and 
rms noises of each image are given in Table [TJ We also 
note that the largest angular scales that the ALMA and 
VLA images should be sensitive to are «25" and 44", 
respectively. 

To measure the integrated flux densities and source 
sizes from each field, we used the task imfit in CASA 
to fit sources within a circular aperture having a radius 
of 15". In Table [2] we list the deconvolved source param¬ 
eters from IMFIT, including the positions of the source 
components, peak brightnesses (S'p), integrated flux den¬ 
sities {Si), deconvolved source sizes (^m x ^m), and cor¬ 
responding brightness temperatures (Tb)- In addition to 
the errors reported by imfit, the uncertainties on the 
peak brightnesses and integrated flux densities include 
the contribution from image rms, as we ll as an absolute 
calib ration uncertainty of 3% at 33 GHz (jPerlev fc Butled 
I2013I1 and 5% at 94.5 GHz (ALMA Cycle 1 Technical 
Handbook), which in fact dominates the uncertainties. 
Clearly for the case of Enuc. 2, which is an elongated 
structure at the end of the bar, a single-component fit is 
not the most appropriate. 


2.5. Line Imaging and Photometry 

We detected the J = 1 ^ 0 lines of HCN and HCO+ 
towards each targeted region in NGC 3627. We did not 
detect H40a for any of the targeted regions, and our fre¬ 
quency coverage missed HNC (J = 1 —> 0) in each case. 
Similar to the continuum imaging, the line data were im¬ 
aged using the task CLEAN in CASA. However, before the 
line images were created, the data were first continuum 
subtracted using the CASA task UVCONTSUB. The rms 
noise of the ALMA channel maps is Ril.2mJybm“^ in a 
10kms“^ channel. 

Moment 0 maps were constructed by integrating the 
spectra at each pixel of the spectral cubes. These are 
overlaid on the continuum maps in Figure Unfor¬ 
tunately, the detection of HCN towards the nucleus is 
cut off by the bandpass at a radial velocity of ulsrk ~ 
800kms“^. The moment 0 map for HCO“'" towards the 
nucleus shows that we are missing higher velocity HCN 
emission to the north. This can also be seen in Figure |3l 
which shows the spectra of both HCN and HCO+ for all 
targeted regions. The spectra were extracted using a cir¬ 
cular aperture centered on the 94.5 GHz source positions 
returned by imfit having a diameter that was 20% larger 
than the corresponding source size major axis convolved 
with the 94.5 GHz beam. We therefore exclude the nu¬ 
clear HCN line emission in this study. For all other re¬ 
gions in Figure m the HCN and HCO+ line profiles are 
fit by a Gaussian. The corresponding velocity dispersions 
are given the each panel, along with uncertainties that 
were estimated using the errors per channel along with 
a Monte Carlo approach. 

Given that the CO(J=I—>-0) image has a much 
coarser angular resolution than our ALMA line images, 
we must first convolve our images to the resolution of 
the CO map before estimating dense molecular gas frac¬ 
tions. This must similarly be done for the continuum (33 
and 94.5 GHz) data for estimating star formation efficien¬ 
cies (see 113.31) . The resolution matching was carried out 
using the task imsmooth in CASA, before re-gridding 
each image to a common pixel scale. For such measure¬ 
ments requiring the convolved data, photometry was car¬ 
ried out by simply taking the peak brightness within a 
circular aperture having a radius of 15" (see Table 13, 
as this provides a way to mitigate any differences in our 
comparisons that may arise from being observed with 
different interferometers having different sensitivities on 
different spatial scales. 

3. RESULTS AND DISCUSSION 
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Table 3 

Photometry Incorporating CO Data 


Source 

Ico 

(KKms-i) 

%CN 

(K Kms-l) 

^HCO+, 

(K Kms-l) 

S 33 GHz 
(mJy bm“^) 

NGC 3627 

375.47± 19.3 

4.54±0.48 

13.80±0.81 

2.10±0.06 

NGC 3627 Enuc. 1 

48.19± 4.4 

1.05±0.23 

0.50±0.20 

1.56±0.05 

NGC 3627 Enuc. 2 

263.03± 13.6 

6.62±0.46 

7.15±0.47 

3.38±0.10 


Note. — Photometry was carried out after convolving both VLA and ALMA 
maps to the 7'.'3 X sythesized beam of the BIMA SONG J = 1 —>• 0 CO data. 



500 600 700 600 900 1000 500 600 700 800 900 1000 

Vlsrk (km s'') 


Figure 3. HCN and HCO”*" spectra (5kms“^ channels) for each 
of the three targeted regions within NGC 3627 extracted using a cir¬ 
cular aperture centered on the 94.5 GHz source positions returned 
by IMFIT. The diameter of each aperture was 20% larger than the 
corresponding source size major axis convolved with the 94.5 GHz 
beam. The full velocity range of HCN emission towards the nucleus 
is cutoff by the bandpass at ulsrk ~ 800kms“^. For all other re¬ 
gions, the spectral lines are fit by Gaussians (dashed lines), and 
the corresponding fitted velocity dispersions and uncertainties are 
given in the upper left corner of each panel. 

In the following section we present our results along 
with a brief discussion about their implications for the 
star-formation activity in the three regions observed in 
NGC 3627. 

3.1. Spectral Indices and Thermal Radio Fractions 

We created a 33-to-94.5 GHz spectral-index map (i.e., 
a, where Si, oc v°‘) of each source after convolving the 
slightly higher resolution 94.5 GHz data to match the 
beam of the 33 GHz data and putting the data on the 


same grid. In the top panels of Figure |4] the spec¬ 
tral indices are overlaid on 94.5 GHz continuum images 
for which pixels not detected at the 3 (t level have been 
clipped. Focusing only on the regions (i.e., pixels) de¬ 
tected at > 3 ct in each map, we measured mean spectral 
indices of = -0.10,-0.15, and —0.27 for the 

nucleus and extranuclear regions 1 and 2, respectively. 
The corresponding standard deviations over these re¬ 
gions are 0.20, 0.18, and 0.15, respectively. As expected, 
these are relatively flat, being < —0.5 for each discrete 
source. We note that the spectral indices in some cases 
become positive in the outskirts of the Hii regions, which 
may indicate a significant contribution of thermal dust 
to the 94.5 GHz emission, although the uncertainties on 
the spectral indices in the Hii region outskirts become 
large. 

If we assume a fixed non-thermal spectral index for 
each source, we can use the measured spectral indices to 
estimate t he fractional contributions from therm al emis¬ 
sion (e.g., iKlein et ^11984 iMurphv et al.ll2012D . While 
most sensitive to the assumption of the non-thermal 
spectral index, this simple thermal decomposition also 
assumes that the free-free emissio n does not become op - 
tically thick at ^ > 30 GHz le.g.. iMurohv et ani2010bll . 
and that there is an insignificant contribu tion of anoma¬ 
lous in icrowave emission at '^33GHz (e.g.. lMuri)hv et all 
l2010al) and thermal dust emission at ^^94.5 GHz. As dis¬ 
cussed below, the few regions where there may be evi¬ 
dence for thermal dust emission are excluded from the 
thermal fraction calculations. We took the non-thermal 
spectral index to be = —0.85, which is the average 
non-thermal spectral index found amon g the 10 star- 
formin g regions studied in NGC 6946 bv iMurohv et al.1 
1 201 Ih and very similar to the average value found by 
iNiklas et al.l (I1997L i.e., = —0.83 with a scatter of 

(Tq,nt = 0.13) globally for a sample of 74 nearby galax¬ 
ies. The 33 GHz thermal fraction contours are overlaid 
on the 33 GHz continuum images in the bottom panels of 
Figure m indicating thermal fractions near unity on the 
peaks of the continuum emission. 

While NGC 3627 is known to host an AGN, the AGN 
does not appear to contri bute significantly to the ra¬ 
dio continuum emission. iFilho et al.l (j2004l) reported 
a 5 GHz flux density of < 0.3 mJy within a 3.6 mas x 
1.4 mas beam using the VLB A. Scaling this to 33 GHz 
assuming a non-thermal spectral index of = —0.85 
results in an estimated AGN contribution to the 33 GHz 
emission of <60 ^Jy , or <3% of the integrated 33 GHz 
emission for the nuclear component reported by imfit 
(see H2.4p . We therefore assume that the nuclear emis¬ 
sion is primarily powered by star formation in the rest of 
the analysis. 
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Figure 4. Top row: Spectral index contours between 33 and 94.5 GHz overlaid on the 94.5 GHz continuum image of the nucleus (left), 
Enuc. 1 (middle), and Enuc. 2 (right). In each case, the spectral indices flatten at the peaks of the continuum emitting regions. Bottom row: 
Estimated thermal fractions at 33 GHz overlaid on the 33 GHz continuum image of the nucleus (left), Enuc. 1 (middle), and Enuc. 2 (right) 
of NGC3627. For both extranuclear star-forming regions, the 33 GHz thermal fractions approach unity on the peaks of the continuum 
emitting regions. 


Ignoring those regions for which the spectral index was 
measured to be flatter than = —0.1 by more than 
Icr, where there may be a contribution from thermal 
dust, we measure corresponding mean 33 GHz thermal 
fractions over the nucleus, Enuc. 1 and 2 of 85, 82 and 
70%, respectively. The corresponding standard devia¬ 
tions measured over these regions are 19, 19 and 20% 
The thermal fraction distribution and ave rage values are 
in agr eement with the results presented bv lMurphv et al.l 
(|2012f) . who reported an average 33 GHz thermal fraction 
of ~76% with a dispersion of 24% for their entire sam¬ 
ple, and even higher values (i.e., >90%, on average) for 
all sources resolved on scales of <500 pc. We note that 
assuming a much flatter non-thermal spectral index has 
little effect on the estimated thermal fraction given the 
flat spectral indices measured between 33 and 94.5 GHz. 
For example, by instead assuming = —0.6, the mean 
33 GHz thermal fractions over the nucleus, Enuc. 1, and 
2 are 80, 76, and 62%, respectively, which are well within 
the quoted standard deviations. 

3.2. Spatial Offsets Between Continuum and Line 
Emission 

One of the most obvious results from an inspection of 
the maps in Figure [2] is the spatial offset between the 
peaks of the continuum and line emission. This is es¬ 
pecially evident for Enuc. 2 (i.e., the bar end) for which 
both dominant HGN and HCO'*' peaks are well to the 
east of the 33 and 94.5 GHz continuum peaks. The con¬ 
tinuum and molecular gas are spatially offset in the plane 


of the sky by f»3"(~130 pc at the distance of NGC 3627). 

The critical density of both the J = 1 —>■ 0 lines of 
HGN and HCO< is ricrit ^ 10®cm“^, and thus probes 
the densest, UV-shielded gas that is likely in the pro¬ 
cess of, or soon will be, forming stars. Assuming that 
the 33 and 94.5 GHz continuum is primarily powered by 
fre e-fre e emission, as suggested by the results presented 
in and thus robustly traces ongoing star formation 
that is <10Myr old, this offset suggests that we have 
reached the spatial resolution at which we can reliably 
separate the fuel stockpiles from the current generation 
o f star form a tion. 

iPan et aH (j201.3[ l similarly report '^100 pc offsets be¬ 
tween peaks of dense gas (HGN J = 1 —>• 0) and star 
formation (3 cm radio continuum) in the starburst ring 
of the barred galaxy NGC 7522. These authors suggest 
that dense-gas formation is promoted by gas cloud col¬ 
lisions at the intersections between the galaxy bar and 
ring orbits, and that offsets between the dense gas peaks 
and those of ongoing star formation occur farther down¬ 
stream from the orbit contact points. Similar conclusions 
have also been found for nuclear bars in nearby starbursts 
te.g.- IMeier eF^[20Ml . Like NGC 7522, N GC 3627 is a 
barre d galaxy with an inner molecular ring (|Regan et al.l 
120021 1 ■ suggesting that similar dynamics may also be driv¬ 
ing the formation of dense gas and its observed spatial 
offsets from the peaks of ongoing star formation activity 
in Enuc. 2. 

Assuming this scenario to be true, we can make a rough 
estimate of the propagation time between the newly 































Microwave Continuum Emission and Dense Gas Tracers in NGC3627 


7 




0.10 

0.08 

0.06 - 

0.04 - 

0.02 - 

0.00 

80 

r- 60 - 


E 


40 - 


o 

u 

X 

b 


20 


2.0 - 


8 1.5 

X 


1.0 


0.5 

0.0 


I 


* 


■ A. 
y*\ 




Nuc. * 
Enuc. 1 O 
Enuc. 2 A 


2 3 

SEE (Gyr'‘) 




-e- 


Figure 5. Top: Dense gas fraction plotted as a function of star 
formation efficiency for each region observed in NGC 3627. The 
quantities on the two axes are both normalized by CO luminosity 
and so some spurious correlation is expected in the presence of 
noise. However, to the contrary we observe an anti-correlation 
between star formation efficiency and dense gas fraction, which is 
opposite the sense expected from only noisy CO measurements. 
Middle: HCO”*" velocity dispersions from the fitted line profiles 
shown in Figure[3]plotted as a function of star formation efficiency 
for each region observed in NGC 3627. Bottom: Ratio of HCN 
to HCO”*" brightnesses plotted against star formation efficiency, 
indicating a trend of increasing line brightness ratio with increasing 
measure of ionizing photons per unit molecular gas mass. Given 
that we are missing higher velocity HCN emission to the north of 
the nucleus, this data point is omitted in the bottom panel. 


formed stars, as traced by the free-free radio continuum 
emission, and the dense gas peaks. The rotation speed of 


the NGG 3627’s disk at the distance to the bar en d (i.e ., 
f»100" Ri4.5kpc) is R:!204kms~^ (|de Blok et al.l 120081) . 
The bar pattern sp eed for NGC 3627 is «55kms“^ kpc“^ 
(jSheth et al.ll200^ . corresponding to f»247kms“^ at this 
distance, yielding a net velocity of between the disk and 
the bar of Ri43kms“^. Thus, the propagation time is 
Ri3 Myr, less than the age of typical OB stars. 

Another striking result from comparing the Enuc. 1 
maps in Eigure [2] is that the HCN and HCO"*" peaks 
are not co-spatial. The peak of the HCO'*' emission is 
located near the peak of the 33 and 94.5 GHz contin¬ 
uum emission. However, the HCN emission peaks in re¬ 
gions that surround the HCO'*' and continuum peaks. 
This may be the result in variations in the (dense) gas 
density, given that the critical density of HCN is nearly 
and order of magnitu de larger than that of HCO'^ (e.g., 
IMeier fc Turned 1201 21) . Stockpiles of ultra dense gas are 
likely still present in the immediate vicinity surrounding 
the newly formed Hii region. Although, excitation effects 
may also be in play (e.g., see 113.3p . 

3.3. Dense Gas Fraetions and Star Formation 
Efficiencies 

In the top panel of Figure [5] we plot the dense molec¬ 
ular gas fraction against the estimated star formation 
efficiency for each region. The dense molecular gas 
fraction is taken to be the mass ratio of dense molec¬ 
ular gas traced by the J = 1 ^ 0 transition of 
HCO'*' to the total molecular gas traced by the J = 
1 —>■ 0 transition of CO. For this calculation we as¬ 
sume HCO'''-to-H 2 and CO-to-H 2 conversion factors of 
Q!unn + ~ 10M(7 i (Kkms~^pc^) ^ le.g.. [Gao fc SolomonI 
[2(TO we assume Q!hco+ ~ Qihcn) and aco ~ 
4.3M(7i (Kkms~^ l^oloiaon fc Vanden Bond 

l200,5l : iBolatto et 120131 ). respectively. We only make 
use of the HCO"*" for deriving the dense molecular gas 
masses as we have data for all regions. The star for¬ 
mation efficiency (SEE), whose inverse is often referred 
to as the gas depletion time, is taken to be the ratio of 
the star formation rate (SFR) to the total molecular gas 
mass traced by the J = 1 0 transition of CO. We con¬ 

verted the 33 GHz spectral luminosity into a star forma¬ 
tion rate using the calibration for p ure thermal emission 
given in (iMurohv et al.l[20llL 1201^ and as sum ed an av¬ 
erage 33 GHz thermal fraction of 79% (see >13.II) . A clear 
trend of decreasing dense gas fraction with increasing 
star formation efficiency is observed. More specifically, 
the nuclear disk of NGC 3627 appears to have the largest 
fraction of its molecular gas in a dense phase, followed by 
the star-forming region at the end of the bar (Enuc. 2), 
and the isolated Hii region (Enuc. 1). However, Enuc. 1 
(the isolated Hii region), appears to be converting its 
molecular gas into stars more efficiently (i.e., ~ a factor 
of 3 more rapidly) than either the star formation region 
at the end of the bar or the nuclear disk. By plotting the 
HCO'*' velocity dispersions from the fitted line profiles in 
Figure[3] against star formation efficiency in the middle 
panel of Figure 0 we find that the star formation effi¬ 
ciency for each region increases with decreasing velocity 
dispersion. Thus, perhaps unsurprisingly, the dynamical 
state of the dense gas appears to have a larger impact 
on the star formation process than the actual fraction of 
dense gas that is available for star formation. 
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Comparing the dense gas fraction and the efficiency o f 
star formation on much larger scales. lUsero et ^ (120151) 
recently found systematic trends in both the dense gas 
fraction and the efficiency with which dense gas forms 
stars. That study included NGC3627, though with a 
beam area two orders of magnitude larger than our study. 
Similar to our results, but using HCN, they found that 
in « 30 disk galaxies, the apparent dense gas fraction 
increases by moving from the outer, low surface density 
parts of galaxies to the inner regions. At the same time, 
they found the apparent efficiency with which dense gas 
forms stars (for them, SFR/HCN) appears to decrease 
as one moves from the lower surface density disk to the 
inner parts of galaxies. They showed that these results 
could be explained by models of turbulent clouds (e .g.. 
IKrumholz fc McKeell200^ iFederrath fc Klessenj[2M^ in 

which the average density and the turbulent Mach num¬ 
ber in a cloud affect both the dense gas fraction and 
the ability of gas at different densities to form stars. 
Similar results have been found considering the appar¬ 
ently low rate at which dense gas forms stars in the cen¬ 
tral part of the Milky Wav fe.g.. fL ongmor e et al.ll^013l : 
iKauffmann et al]l2013l : IRathborne et al.ll2014fl \^th the 
improved resolution offered by ALMA, we are able to 
measure a velocity dispersion that may be more directly 
related to the turbulent velocity dispersion (at least out¬ 
side the central region). The observation of variable 
dense gas efficiency (SFR/HCO+) and suppressed star 
formation in the central, high dispersion part of the 
galaxy appears to agree at least qualitatively with these 
results. 

Since we have relied on only the J = 1 ^ 0 transition 
of HCO+ to derive the dense gas mass, and have assumed 
similar conversion factors for HCN and HCO+ to dense 
H 2 , it is illustrative to see how the HCN-to-HCO“*' line 
brightness ratio varies among regions. This is shown in 
the bottom panel of Figure [SI and plotted against star 
formation efficiency. The HCN-to-HCO’*' line brightness 
changes by a factor of ^2 between both extranuclear star¬ 
forming regions over a similar change in star formation 
efficiency. This trend may be due to gas excitation ef¬ 
fects given that the abscissa is a measure of the number 
of ionizing photons (free electrons) per unit molecular ga s 
mass. For instance, as discussed in lPapadopouloj (1200711 . 
the HCO+ abundance is known to be sensitive to the ion¬ 
ization degree of molecular gas, which can significantly 
reduce the HCO’*' abundance in star-forming and highly 
turbulent molecular gas, while HCN remains abundant. 
Given the observed trend, this appears to be a plausible 
explanation. However, there are many alternative expla¬ 
nations for changes in the ratio of HCN-to-HCO+ line 
brig htness ratio including variations in gas density (see, 
e.g., iMeier fc Tnrnei]l2012fl . Given the large error bar on 
the line-brightness ratio for Enuc. 1, further speculation 
here is unwarranted. Regardless, it is worth pointing out 
that even if one ignores the large error bar, the trend 
observed in the top panel of Figure [5] would persist even 
if we instead used the HCN-derived dense gas mass. 

4. CONCLUSIONS 

In this paper we have combined ALMA/Band-3 line 
and continuum images with VLA/Ka-band images to 
characterize the star formation activity on R::!l00pc scales 
around three distinct regions within the nearby galaxy 


NGC3627; i.e., a nuclear disk hosting an AGN, an iso¬ 
lated Hii region (Enuc. 1), and a star-forming complex 
at the end of a bar (Enuc. 2). Our conclusions can be 
summarized as follows: 

• The thermal fraction at 33 GHz is nearly unity at 
the peaks of the Hii regions as mapped on «100pc 
scales, with an average value of «76% for both ex¬ 
tranuclear star-forming regions. The mean thermal 
fraction at 33 GHz among the three regions stud¬ 
ied is 79% with a dispersion of 19%. We addi¬ 
tionally found the radio spectral index to become 
both increasingly negative and positive away from 
the peaks of the Hii regions, indicating an increase 
of extended non-thermal emission from diffusing 
cosmic-ray electrons and the possible presence of 
cold dust, respectively. 

• The peaks of the J = 1 —>■ 0 HGN and HGO+ line 
emission are spatially offset from the peaks of the 
33 and 94.5 GHz continuum emission. For Enuc. 2, 
the continuum and gas are spatially offset in the 
plane of the sky by R::i3"(~130pc at the distance of 
NGG3627). Assuming that the 33 and 94.5 GHz 
continuum is primarily powered by free-free emis¬ 
sion, and that the J = 1 —>■ 0 HGN and HCO'*' line 
emission is tracing dense, UV-shielded gas, this in¬ 
dicates that our data reach an angular resolution 
at which one can spatially distinguish sites of cur¬ 
rent star formation from the fuel stockpiles for the 
next-generation of star formation. 

• Combining our ALMA and VLA observations with 
archival BIMA CO (J = 1 —0) data, we calculate 
dense gas fractions and star formation efficiencies 
for each region, finding that the dense gas frac¬ 
tion decreases with increasing star formation effi¬ 
ciency. This suggests that an increase in the dense 
gas content of star-forming regions does not reflect 
an increased efficiency for which parts of galaxies 
can turn molecular gas into stars among the three, 
diverse regions studied here. Specifically, Enuc. 1 
(an isolated Hii region) appears to be converting 
molecular gas into stars more efficiently than both 
the star-forming region at the end of the bar or the 
nuclear disk. We additionally find that the velocity 
dispersion of the dense gas in each region decreases 
with increasing star formation efficiency, indicating 
that the dynamical state of the dense gas, rather 
than its abundance, plays a larger role affecting on¬ 
going star formation activity. 

• The J = 1 —0 HCN/HCO’*' brightness ratio and 
star formation efficiency both vary by a factor of w2 
between the extranuclear Hii regions. We speculate 
that this may be due to an increase in the ionization 
degree of molecular gas for the star-forming region 
having a higher star-formation efficiency, which in 
turn reduces the HCO+ abundance while HCN re¬ 
mains abundant. However, alternative explana¬ 
tions (e.g., variations in the gas density) may also 
be plausible. 
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